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Abstract 

We study the distribution of the n-th energy level for two different one-dimensional 
random potentials. This distribution is shown to be related to the distribution of the distance 
between two consecutive nodes of the wave function. 

We first consider the case of a white noise potential and study the distributions of energy 
level both in the positive and the negative part of the spectrum. It is demonstrated that, 
in the limit of a large system (L — > oo), the distribution of the n-th energy level is given 
by a scaling law which is shown to be related to the extreme value statistics of a set of 
independent variables. 

In the second part we consider the case of a supersymmetric random Hamiltonian (po- 
tential V{x) = (j>{x)'^ + 4)'{x)). We study first the case of <t>{x) being a white noise with 
zero mean. It is in particular shown that the ground state energy, which behaves on average 
like exp —L^^'^ in agreement with previous work, is not a self averaging quantity in the limit 
L — > oo as is seen in the case of diagonal disorder. Then we consider the case when (j){x) has 
a non zero mean value. 

1 Introduction 

One-dimensional disordered systems have been studied in great detail in the past and are still a 
subject of interest. Examples include recent work on the failing of the single parameter scaling 
for localization near the band edge or at strong disorder the statistical properties of the 
Wigner time delay, studied in several works |3|, ^, |6|, as well as the correlations of the 
time delay for different energies Q. The distribution of another transport time has also been 
investigated in As a last example let us mention that the ac-conductivity was re-examined 
ini- 

In this article we are interested in the following one-dimensional Hamiltonian: 

^ = -^ + ^(^)' (1) 

where V{x) is a random potential (we choose units such that h = 2m = 1 for simplicity). The 
spectral properties of such Hamiltonians were studied in different works for various kinds of 
disorder. The case of (5-scatterers of random weights and/or at random positions was examined 
by several authors |1C, 11, |l^ (see also Here we consider two kinds of random potentials. 



(A) In the first case the potential is a white noise, i.e. V{x) is distributed with the Gaussian 
weight: 

VV{x) P[V{x)] = VV{x) exp - — / dx V{xf . (2) 

2a J 
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This implies that {V{x)) = and {V{x)V{x')) = a 6{x — x'), all other cumulants being zero. The 



spectral properties of this model were studied by Halperin |14]; later Berezinskii developed a 
diagrammatic method adapted to this particular model used to study various quantities in a 
number of works. We can also mention a recent study of the density of states and level statistics 
of this model using soliton calculus This model is equivalent to the 5-scatterer model in a 
certain regime: when the density of impurities is very large compared to their weights. Moreover 
it was shown in |]l^ that results for the Gaussian disorder are very generic at high energy: they 
reproduce what is expected for any random potential provided that the correlation length of this 
disordered potential (which is zero for the Gaussian disorder) and the de Broglie wavelength are 
the smallest length scales of the problem, whatever their relative magnitude is. 

(B) The second disordered potential we will consider belongs to another class of random 
potentials, which has attracted the attention of several authors: the case of off-diagonal disorder. 
The random potential reads 

V{x) = 4>{xf + (P'{x) (3) 
where ^(2;) is random. The Hamiltonian, that may be factorized as H = Q'^Q where Q = 



—dx + 4>{x), describes one-dimensional supersymmetric quantum mechanics ]18|]. A discrete 
version of this model would be a ID tight-binding Hamiltonian with random hopping whereas 
the discrete version of model A is a ID tight-binding Hamiltonian with random site energies. 
The interest for Anderson model with off-diagonal disorder was recently renewed due to its 



connection with disordered spin chain models |1£, It is also worth mentioning the 



relation of the supersymmetric Hamiltonian with the problem of classical diffusion in random 



media (see for example [22] and references therein); many properties were obtained in this 
context very recently using a real space renormalization group method |^3|. In the case of 
off-diagonal disorder, the density of states presents different kinds of behaviour compared to 
diagonal disorder, like Dyson singularities at band edge [U] (if {4>{x)) = 0), or power law 
singularities (if {(l){x)) / 0) |2|, M, pi p|, |0|, P, ||1 . Very recently the density of states 



for coupled chains with off-diagonal disorder was studied |33| showing interesting effect of the 



parity of the number of chains. For a complete review on one-dimensional disorder systems. 



the interested reader is refered to [13, In the following we will consider the case where 
the random function involved in the potential is distributed according to a Gaussian weight: 
V(^{x) P[cP{x)] = V<P{x) exp / dx [(t){x) - ^^gf . 

The purpose of this article is to study the distribution of the energy En of the n-th excited 
state of the Hamiltonian (0) (n + 1-th energy level) considered on a finite interval of length L 
with the wave function '^{x) satisfying Dirichlet boundary conditions: (/^(O) = <^(L) = 0. Let us 
denote this distribution: 

Wn{E) = {5{E - En)) , (4) 

where the brackets (• • •) mean average over different configurations of the disordered potential 
V{x) (with respect to measure (^) for model A). These distributions are related to the average 
density of states per unit length by the relation: 

^ 00 

p{E) = -Y,Wn{E) . (5) 

n=0 

This problem was studied by Grenkova et al. who derived these distributions for the 
(5-impurity model of Frisch and Lloyd when the weights of impurities are very large compared 
to the average density of impurity and the energy. We stress that this is not the limit where 
this model is equivalent to the model with the white noise potential considered here. 

The problem was later addressed by McKean [36[ who derived the distribution for the ground 
state energy (n = 0) for < 0; he considered different boundary conditions for the model 
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originally studied by Halperin for a potential being a white noise. We will provide in section |3| 
a generalization of the result of McKean by giving the distribution for all eigenvalues in both 
regions of the spectrum E < and E > 0, provided \E\ ^ cj^/^. We will demonstrate that the 
distribution is given by a scaling law in the limit L ^ oo and will show how the parameters 
scale with the sample size L. This scaling law is similar to the distribution of the extreme value 
of a set of statistically independent random variables |37]. This is in agreement with the fact 
that the eigenvalues are not expected to present level repulsion in the limit for which the states 
are localized, as demonstrated by Molcanov [3^ . 

After having considered the case A of diagonal disorder we will study in section ^ the problem 
for off-diagonal disorder B at band edge and at high energy as well. We will consider first the 
case fj, = for which we will study the ground state energy distribution. The analysis for model 
A cannot be applied in this case and we will need to use specific approximations. Note already 
that our result for the distribution of the ground state energy gives a mean value in agreement 
with the prediction of Monthus et al. |39| who showed that the averaged ground state energy 
behaves like if the system size is very large, by finding a lower bound and an upper bound. 

Moreover, our result shows that the ground state energy is not a self averaging quantity as for 
model A at low energy, and its distribution presents a large tail. We also give the distribution 
Wn{E) in the high energy limit. We finally study the distributions for /x 7^ in the low energy 
limit. 



2 Idea of the method 

In this section we give the main ideas of the method we use to derive the eigenvalue distribution. 
We concentrate on model A since the ideas are the same for model B apart from small subtleties 
which we will discuss later. 

Let us consider iIj{x; E), the solution of the Schrodinger equation Hip{x; E) = Eip^x; E) with 
the boundary condition ip{0;E) = 0. The boundary condition 'ip{L;E) = is fulfilled whenever 
the energy E coincides with an eigenvalue En of the Hamiltonian. In this case, the wave function 

^nix) = tp{x; En) I /q dx' ^(x'\ Enf' has n nodes in the interval ]0, L[, and two nodes at the 
boundaries. Let us denote by the n + 1 lengths between the nodes. We consider the Ricatti 
variable 

z{x-E) = ^^\^\^[x-E)\ , (6) 
which obeys the following equation: 

^z = -E- ^V{x) , (7) 

for an initial condition z(^\ E) = +00. This equation may be viewed as a Langevin equation for 
a particle located at z submitted to a force deriving from the unbounded potential 

U{z) = Ez + j (8) 

and to a random "force" V{x) (white noise). Each node of the wave function corresponds to 
\z{x)\ = 00. At "time" x = the "particle" starts from z{0) = +00 and eventually ends 
at z{ii — 0"*") = —00 after a "time" ii. Just after the first node it then starts again from 
z{£i + 0^) = +00, due to the continuity of the wave function. It follows from this picture that 
the distance im between two consecutive nodes may be viewed as the "time" needed by the 
particle to go through the interval ] — 00, +cxd[ (the "particle" is emitted from z = +00 at initial 
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"time" and absorbed when is reaches z = — oo). Fohowing appendix^ we introduce the n-th 
moment Cn{z) of the "time" C the particle takes to reach — oo starting from z: 

Cn{z) = (£" I z(0) = z; z{C) = -oo) . (9) 



According to ( |180D or ( 188 ) these moments satisfy the following recursion relations: 



z „ r+oo 



= ^/ Az'e^uiAl dz"e-^^(^") (10) 



U{z) = 1. (11) 

The boundary condition at a = +oo is chosen to be a reflecting one. This choice, which simplifies 
calculations, is not important since the particle can never go back to this edge of the interval 
and necessarily ends at h = — oo. 

The moments of the lengths between consecutive nodes of the wave function are: 

(r) = £„,(+oo) . (12) 

It is worth mentioning that the n + \ lengths are statistically independent because each time the 
variable z reaches — oo, it loses the memory of its earlier history since it is brought back to the 
same initial condition and V{x) is (5-correlated. This remark is a crucial point for the derivation 

0iWn{E). 

However the fact that V{x) is (5-correlated is not essential to ensure the statistical indepen- 
dence of the ^m's. Indeed, imagine that correlations of V{x) are short range, on a scale Xc- 
Equation (^) shows that when z starts from +oo at initial "time" , it needs a "time" Ax to reach 
the region in z-space where U{z) is of order or smaller than a and where the presence of the 
random "force" V{x) matters for the evolution of z. The "time" Ax, during which the dynamics 
of z is governed only by the deterministic "force" —U'{z), can be defined as Ax = e+z'^ 
where U{zq) = a; we have Ax ~ l//c for \E\ S> cr^/^, and Ax ~ a~^^^ for \E\ <C <t^/^. If z follows 
the deterministic evolution from +oo during a "time" Ax sufficient for V{x) to decorrelate, then 
the £m's are statistically independent; this occurs when Xc is much smaller than the smallest 
length scale among 1/k and cj"^/^. 

Let us finally write the equation satisfied by the generating function of the moments of the 
traversal "time" C 

h{a, z) = (e-°^ I z(0) = z; z{C) = -oo) . (13) 



According to ( |183D it obeys: 

Gz h{a, z) = a h{a, z) , (14) 

where 

= -U'{z)dz + (15) 

is the generator of the backward Fokker Planck equation (BFPE) associated with the stochastic 
differential equation —U'{z) = —dzU{z) is the force deriving from the potential. The 

generating function satisfies the boundary conditions: 

dzh{a,z)\z=+oo = 0, (16) 
/i(a,-oo) = 1. (17) 

Coming back to the initial problem, our goal is to compute the probability for the energy 
of the n-th excited state to be E. This occurs if the sum of the n + 1 distances between the 
nodes is equal to the length of the system: L = Yl^=i ^m- As it was stated above the £m's 
are independent variables and Prob L = is given in terms of the distribution of the 

variables im- For the different cases we will analyze throughout this article, we will initially 
examine the distribution P{t)^ enabling us to find Wn{E). 
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3 Diagonal disorder: white noise potential 



3.1 Distribution of the distance between consecutive nodes of the wave func- 
tion 

We study the distribution P{i) of the distance i between two consecutive nodes of the wave 
function. Let us first note that the average integrated density of states per unit length N{E) = 
dE' p{E') gives the number of states below energy E per unit length, which is also the 
number of nodes of the wave function of energy E per unit length, or in other words the inverse 
of the average distance between two consecutive nodes: 



{£) =Ci{+oo) = NiE)-^ . 

The average integrated density of states is given by calculating (|l^) for n 
possible to perform integration over z" |14, |T^ ]: 



N{E) 



(a/2)V3 



'IT 

7r2 



dy 



^(M u 

V 12 ^ 



Ai^ 



-E 



:a/2)2/3 



((t/2)2/3 



+ Bi^ 



y) 



:ct/2)2/3 



(18) 

1 in which it is 

(19) 
(20) 



Ai(x) and Bi(a;) are Airy functions. This result, given by Halperin |l^, was first mentioned 



in |11] as an approximation for the integrated density of states for a potential consisting of 
randomly dropped (5-scatterers in the limit of high density of scatterers. 
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Figure 1: Density of states |14] per unit length for cr = 1, given by (19). Dashed: free density 
of states. 

We now examine the limit of high energy \E\/a'^^^ — > cxd in the negative and in the positive 
part of the spectrum. 



3.1.1 Negative part of the spectrum E = —k^: trapping of the Ricatti variable 

li E < the potential U{z) possesses a local minimum at z = \/—E = k that may trap the 
Ricatti variable (see figure §). In the limit \E\/a'^/^ ^ 1, the well is very deep (or the diffusion 
very small) and when traveling from +oo to — oo, the Ricatti variable spends most of the "time" 
in the well. Then we expect that the average "time" {£) is given by the Arrhenius law and 
that its distribution is a Poisson distribution, as demonstrated in appendix Expanding the 
potential U {z) in the neighbourhood of its two local extrema 

U{z) ~ - +k{z-kf (21) 
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Figure 2: The potential that traps the Ricatti variable z for a = \ and E = —4. 



U{z) 



2k^ 

-k 3 



(22) 



equation ( |193| ) gives 



{l) = Ci{+oo) = N{E)-^ ~ 



vr 



3/2 



(23) 



This argument was first used by Jona-Lasinio |40| to find the exponential factor of the integrated 
density of states. We stress here that it also gives the correct pre-factor, that may be checked 
by extracting the limiting behaviour of (19). 

As we have shown in appendix ^ the distribution of the time spent in the well is a Poisson 
law: 

P{e) = N{E) exp -^ N{E) . (24) 



This equation will be the starting point to find the distribution Wn{E) in section 3.2. 



3.1.2 Positive part of the spectrum E = +fc^: small disorder expansion 

In the positive part of the spectrum, the potential in which the Ricatti variable evolves has no 
local extremum (see figure ^) and therefore the length (£) does not anymore follow an Arrhenius 
law. 




Figure 3: Potential for 

In the absence of diffusion (cr = 0) the "time" 
to — oo is i = n/k, which may be found either by 



a = 1 and E = +4. 

needed by the Ricatti variable to go from +cxd 
integrating (0) in the absence of the potential 
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V{x) or by taking the limit cr ^ in (18,19[). We can expect that, for sufficiently weak disorder, 
the length is weakly fluctuating around its mean value. We are now going to show that the 
distribution of £ is a narrow Gaussian distribution in this limit. For this purpose we will analyze 
the moments = Cn{+oo), which is more conveniently achieved by studying a generating 
function. Instead of the generating function of the moments ([l^), it is more advantageous to 
consider the generating function of the cumulants: 

w{a, z) = In h{a, z) . (25) 

It follows from (|l|) and (|l^ that 

- (k^ + z'^)d,w + I [d^w + id.wf] = a . (26) 
The solution may be obtained through a perturbative expansion in the parameter a. We write 

(27) 

-oo) = 0. The 0-th order of 



w = w^^) +w^'^^ + w^^) H 



where w^"'^ = ©(cr"). The boundary condition for w is: w{a,z 
m gives 



w^°\a,z) 



-a 



,1 a / Z TT 

,^"FTi^ = -ir^*^"^ + 2 



(28) 



The integral is the "time" needed to go from z to — oo in the absence of diffusion as it is clear 
from (j^. The n-th order of (|2^) gives: 



a 



w^'^\a,z) = - 



dz' 



fc2 + z'2 



n-1 



92,u;("-i)(a,z') + ^ d,m^"'\a,z')d,,w^''-^-'^\a,z') 



m=0 



(29) 



In particular we have: 

w^^\a, z) 



a 

2^ 



a 



2fc (1 + z2/A;2)2 



a\2 Ax 
+ 1^ 



kJ 



z/k 



dx- 



(1 + X2)3 

3 f^/^ 2dx 



(30) 



(l + a;2)5 \kJ 4(l + z2/A;2)4 \k 



a 



(1+X2)5 



(31) 

Equation ( |30| ) gives the dominant contribution to the second cumulant {{(^)) of the length. 
Explicit computation of terms of w of higher order in a is not required and we only need to 
know their behaviour with a. Using a recursion argument it is possible to show that lo^"^ is a 
polynomial of degree n + 1 in a: 



n+l 



W 



in), 



a,z) = a" ^ a"'c„(z) . 



(32) 



m=l 



We now extract from these expressions the information about the cumulants of 



w(q, +00) = 



-a 



m=l 



-an) ■ 



Using (|2|,|0D we see that: 

w^a, +00) 



(a, +00) + w(^) (a, +00) + O(a^) 
vr /7r\2 3a 2^ 



'"A; ^ 2! \kJ SvrfcS 



(33) 

(34) 
(35) 
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then 

it) = l + 0{a') (36) 

The n-th cumulant is given by the term proportional to in w{a, oo). Equation (^) shows 
that the term of lowest order in a containing a" is w^^~^\ then 

{{r)) = 0((7"-i) . (38) 

The n-th order cumulant for n > 2 characterizes fluctuations that are negligible, in the high 
energy limit, compared to the fluctuations described by the second cumulant: 



((£2))"/2 



--1 
2 ^ 



— > 0. (39) 



Since the second cumulant is dominating it follows that the distribution of i is Gaussian in this 
limit: 



For the positive part of the spectrum, ( |40| ) shows that the fluctuations of I are small compared 
to its average value, in contrast to what happens in the negative part of the spectrum where 
(pi) shows that the fluctuations of I are of the same order as the average value. 



The structure of the w^ave function 

Since we are dealing with one-dimensional disordered system, the wave functions are expected to 
be localized, i.e. decreasing with an exponential damping on a length scale being by definition 
the localization length A. Let us recall that A ~ — for E — > +oo and A ~ / ^ for E — > — oo 



(see 1 17] for example). The results we have derived for the distribution P{1) show that in the 
limit E +00, the consecutive nodes of the wave function are separated by weakly fluctuating 
distances of order vr/fc, which is much smaller than A. In contrast, in the limit E — > — oo the 
distance i is distributed according to a Poisson law ( p^ ) which means that it is probable to find 
two consecutive nodes as close as possible. However the typical scale of the length is much larger 
than the localization length: i ~ N{Ey^ ~ Avrexp ^{-E)^/'^ > A. 

3.2 Distribution of individual energy level 

We now derive the distribution Wn{E). As we have mentioned above, E coincides with the 
eigenvalue En if the sum of the n + 1 lengths between nodes is equal to the length L. 

3.2.1 Low energy E = —k^ 

We follow here the idea McKean used to find the ground state energy distribution [^]. The 
probability that the energy E is between two consecutive energies is: 

Prob [En^i <E<En]= Prob [^i + •••+£„< L < + ••• + 4 + (41) 

= / dh P{ii) ■■■ d£n+l P{ln+l) e{ll + ---+ In+l - L) {L - h £„) ,(42) 

Jo Jo 
where 0{x) is the Heaviside function. Using ( p^ we get: 

Prob [En-, <E< En] = ^-Mjpl e-^^(^) . (43) 
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Introducing the joint distribution for the eigenvalues 

T^o,l,..,n,..(^0,^l, ••• •••) = ( n ^^^^ - ) (^^) 



j AX^dXi - ■ -dXn - ■ ■Wox-,nAXo,Xu - ■ ■ ,Xn, - ■ ■) e{E-Xn-i) d{Xn-E) (45) 



\m=0 

and differentiating (pT|) with respect to E gives 

^ Prob <E<En] 

_ _d_ 

~ dE 

Differentiation of the integrand of (|45| ) gives two terms. In each term the Heaviside function 
0{Xn — Xn-i) does not play any role since the joint probability for the different energies is 
proportional to the following product of Heaviside functions: 

VFo,i,...,n,..(^o, Xi, • • • , x„, • • •) a e{Xi - Xo) e{X2 - Xi) • • • e{Xn - x^-i) ■■■ . (46) 

Then we find: 

^ Prob [En-l <E<En]= Wn-l{E) - Wn{E) . (47) 
For n = we recover the result of McKean |3G] for the ground state energy distribution: 

dE'Wo{E')=e-^^^^y (48) 



rrx 

Je 



Using (43) and (|47|), it is now easy to show that: 

n\ 

This result has a clear meaning since Lp{E) gives the probability to find an energy at E^ the 
factor ^e~^ "compelling" the number of states below E, x = LN(E), to be close to n. 
Let us remark that the relation (^) is satisfied. 

We now analyze this result in more detail and demonstrate that the distribution ( ^9|) has a 
scaling form in the limit L — > oo. We use the approximated expression ( p3D of the integrated 
density of states per unit length to write the distribution Wn{E) as: 

n! 

with 



g{E) = {n + l f^~^> + 2(n + l)L^/^e-i(-^)''' - H-E) , (51) 

3 2 

where we have introduced: 

i=4^. (52) 

We set (7 = 1 for simplicity since it will be easy to recover its dependence at the end (the 
dimension of a is: [a] = length"'^ = energy^/^). 

We first look for the typical value of the n-th excited state energy E^'^ . It is the solution of 
the equation g'{E) = 0. The variable Y = ^{-E^^^f/'^ is the solution of 

y _ n+2 y/3 
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In the limit L ^ oo we find Y = 31n(L) + ln(31n(L)) + O 

2 ln(3 In L) 



ln(3 ln(L)) 
ln(L) 



that is: 



ETiL) 



1 + 



9 InL 



+ 



ln2(31nL) 
In^L 



(54) 



where we keep the first correction to the dominant term in lia^^^ L because it is larger than the 
difference between the typical values of two consecutive levels ([6^). This behaviour was given 
as a good approximation of Eo by McKean ||3^. Note that the n-dependence enters only in L. 
We will need the maximum value of the distribution, which is: 



2(n + l)"+ie-("+i) / ~\i/3 



n! 



3 InL 



1 + 



ln(3 In L) 
InL 



(55) 



We now study the derivatives of g{E) at E = E^'^. By definition the first derivative vanishes: 

'ln(31nL)\ 



O 



After a little of algebra, we find for the m-th derivative: 



(56) 



^{m)(-^^yp) = 2"^(n + l) (3 InL 



_ \ m/3 



1 + 



ln(3 In L) 
InL 



m>2. 



(57) 



The second derivative defines the scale of the fluctuations of Er, : 



SEn 



a 



2/3 



.^-1/3 
3 InL 



2^/^^TT 



1 + 



ln(3 In L) 
InL 



(58) 



Note that the width of Wn{E) decreases as L — > 00. The relative fluctuations tends to zero, 
that is En is self averaging in this limit. 

Moreover it is possible to introduce the function UniX): 



WniE) 



1 



UJr. 



with 



MX) 



5En 
Wn{E^r[') 



E - ET 

5Er, 



(59) 



: exp 



^-^ ml 



(60) 



m=l 



where am = g^'^^E'^'^) /[g'^^'^ {E]^'^)]'^/^ . Using (||,|3) it is easy to see that oi = and a„ 
(n + l)i-"^/2 in the limit L ^ 00. Then we flnd for the scaling function: 



UJniX) 



(n + 1) 



exp (Vn + IX - (n + 1) 



(61) 



A similar result was obtained by Grenkova et al. |35] for the (^-impurity model of Prisch and 
Lloyd, in the limit of low impurity density which has no equivalent in the model we consider 
here. These authors found the same scaling distribution (|6l| ) whereas the scaling (|5^j58|) is 
different (it is of course model-dependent). 

It is interesting to note that (61) is related to extreme value statistics. If we consider a set 
{xk} of AA ^ 00 statistically independent variables distributed according to the same law and 



10 



order them: xo < xi < X2 < ■ ■ ■, then the distribution of Xn has the form of (^) up to a rescahng. 
coo{X) is the distribution of the most negative, etc. This problem was studied by E. Gumbel 
in 1935. For more details about extreme value statistics, see [41, 37] (see appendix]^). The 
distribution ( |6T1) is the extreme value distribution for any distribution of the so called exponential 
type (unlimited domain and finite moments). Moreover the form of the scaling ( p^j58| ) allows 
us to find the tail of the distribution p{E) of one of those variables: p{E) ~ exp— | 
This tail is not surprising since it is precisely the behaviour of the density of states (see also 
the discussion in the conclusion). The remark that the eigenvalues are distributed as a set of 
statistically independent variables is in agreement with the expected absence of level repulsion 
for one-dimensional disordered systems in the localized regime ||3^, ^ . Equation (^) allows us 
to compute the generating function 



Gn{k)= I ^ dXu^{X) e^^ = ^^\ n + 1)"^^ ■ (62) 

J-oo r(n + l) 

The expansion of the logarithm of the generating function in powers of k gives the cumulants of 
the distribution (|6l|): 

{X) = ^/^TTT [V'(f^ + 1) - ln(n + 1)] (63) 
((X™)) = (n + l)'™/2^(™-i)(n + l) form ^2, (64) 

where il^iz) = ^lnr(z) is the digamma function. It is possible to show that (^Tj) converges 
to a Gaussian distribution in the large n limit. Using the Stirling formula we find ln^„(A;) ~ 



1 x"^ 

lim uJn{X) = -^=e~~ . (65) 



+ \ ' 

lim LOr,{X) = — 

It is also possible to compute the mean level spacing (A^) between the n-th excited state 
and the n + 1-th one. Equation (^) shows that: 

Kl. - ET = cr'^'l In ^ (3 In l) , (66) 

where L is still given by (^) for the n-th excited state (and not the length associated with the 
n + 1-th one). Since the mean energy is {En) = E^'^ + 6En {X), we eventually find: 

(A„) = {En+i - En) = ^^^^ (3 In l) . (67) 

This result shows that the average distance between two consecutive levels is of the same order 
as the fluctuations ( |58[ ) of the position of those levels. In other words, the distributions Wn{E) 
are overlapping functions, as represented in figure |^. Now, if we consider that L is fixed, we can 
see from ( |58|) and ( |67|) that, apart from the unimportant n-dependence of L in the logarithm, the 
mean level spacing decreases like (A„) oc 1/n whereas the fluctuations decreases like 5 En oc l/\/n 
for large n. As n increases for fixed L, the consecutive distributions Wn{E) become more and 
more overlapped, this is also suggested in figure ^ 

3.2.2 High energy E = 

In order to find the distribution Wn{E) we use the same idea as before: the probability that 
E = En is related to the probability that the sum of the n + 1 lengths between the nodes of 
the wavG function is L. Let us introduce A — ^^rn=i ^m* 

We have Wn{E)dE = P„(A)dA, where 
Vn{X) is the distribution function for A. Since the ^m's are statistically independent Vn{^) is 
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Figure 4: The distributions for the first five energy levels, a = 1 and L = 10 



easily found from or (|^) (note that VqW = = A))- If = En, the sum of the lengths 
between nodes coincides with L: A = L. In this case we have KN{E) = n + 1. Differentiating 
this latter equation we get N {E)dK + K p{E)dE = 0, that is = -j^y- It follows that: 

Wn{E) = ^^Vn{A = L). (68) 



Using (|24D, this equation allows us to get straightforwardly (|49|). In the positive part of the 
spectrum, since the distribution ( ^0|) is Gaussian, the distribution 'Pn(A) is also Gaussian and 
we get: 

Lp{E) {LN{E)-n-lf 



with {{(!. )) being given by (|37|). 

We now go to the variable k = ^/E for clarity. We write: 

Wn{E)d.E = dk^=== exp ; ../^^ , 70 

with s{k) = (n + 1)^ ^^^^ . The maximum value of the distribution corresponds to A; = A;^ = 
(n + 1)^, which is the de Broglie wavelength in the absence of the disordered potential. We 
remember that we are dealing with a high energy limit E ~ ^^^p^ > (t2/3, which may be 
conveniently written as: 

(71) 

where we have introduced the localization length A„ associated with an energy E^ = (kn)'^ = 
(n + 1)2^. We recall that A = ^ for ^ > ct^/^ |l^. The condition is fulfilled either for a 



delocalized regime L <C A, or for a localized regime L ^ A provided that n is sufficiently large. 
Note that there is no restriction on L for the derivation of Wn (E) . 

Due to the condition of validity ([fll), it is possible to neglect the dependence of s{k) on k 
and replace s{k) by s(A;q); then the distribution Wn{E) is a Gaussian distribution 

Wn{E) = ^==e (72) 

of mean: ^ 

{En) = in + lf^ (73) 
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and width 




5El 

{EuV 



According to ( fflD the relative fluctuations are necessarily small: 
interesting is to compare these fluctuations to the mean level spacing (A, 
2(n + Then 



(74) 

^^-(in^i; « 1- More 

= {En+l - En) ~ 



(An) 



3_L_ 

7r2 An 



(75) 



This condition tells us that the fluctuations 5En of the position of En are larger than the mean 
level spacing (A„) if we are in a localized regime; the distributions Wn{E) are overlapping 



functions. This agrees with the fact that no level repulsion is expected in this regime [38, 35 
where the level spacing is believed to be distributed according to a Poisson law [38|. In the 
delocalized regime L <C A, the fact that the distributions Wn{E) are non-overlapping indicate 
level repulsion. 

It is worth mentioning that (73,^) may be found by a simpler, although less systematic, 
perturbative argument: the energy of the n + 1-th level is, up to first order pertubation theory 

En — En + (V'n 1^(^)1 V'n )i with ^J^(x) = Y^sinfcJJx. Then it is straightforward to see that 

{En) is given by (^ and that 6E^ = ([/^f dx F(xX(x)2]2) = a f^dxtl^Kx)^ leads to (||). 

Let us end the section with a remark concerning the ground state energy. We have noticed 
that the ground state energy is a self averaging quantity {{Eq) S> SEq). Its behaviour with L is 
Eq ^ — In^/^ L for L ^ oo, then it is vanishing Eq ^ for a size L ~ cr^^/^ and behaving like 
Eq ~ 1/L2 for L ^ 0. 



4 Off-diagonal disorder: supersymmetric random Hamiltonian 



The Hamiltonian we consider in this section is the following supersymmetric Hamiltonian: 

d2 



Hs 



+ (p{x)^ + (/)' (x) . 



(76) 



The spectral and localization properties of this Hamiltonian were studied for various kinds of 



disorder. The case of white noise was analyzed in |25, 26, 22 1. The case of a random telegraph 
process was also studied in detail |27, Note that the Hamiltonian ( [7^ ) is the square of a Dirac 
Hamiltonian with a random mass. The spectral properties of ID random Dirac Hamiltonians 



were studied in a more general situation by Bocquet [31, 32[. His analysis was extended very 
recently to study the distribution function of the local density of states [42[; these authors used 
replica trick and supersymmetry and obtained for the supersymmetric Hamiltonian in a high 



energy limit a result similar to the one derived by Altshuler and Prigodin [43[ for model A using 
the Berezinskii technique. 

As it was mentioned in the introduction, we focus on the situation where (i>(x) is a white 



noise: 



V(t)ix) P[Mx)] = VMx) exp — 3- / dx Wx) 

25 J 

We recall that the integrated density of states is in this case [^, \ 

N{E) = ^ ^ ^ , 

vr2j2(/B/5) + Ar2(/B/g) ' 

where J^(z) and Ny{z) are the Bessel functions of first and second kind, respectively. 



(77) 



(78) 
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We now consider the situation where = and will discuss the case / in a last section. 



In general, the Hamiltonian (76) may possess a zero mode, however, for the problem of interest 
here, it can not satisfy the Dirichlet boundary conditions; we say that the supersymmetry is 
broken. Due to the super symmetric structure of the Hamiltonian the spectrum of Hs is positive. 
It is possible to rewrite the Schrodinger equation Hs^{x) = k'^ip{x) as two coupled first order 
differential equations: 

Qhix) = k^ix) (79) 
Qifix) = kx{x), (80) 

where Q = —dx + </>(x) and = dx + 4>{x). As in the previous section, the first step is to study 
the statistics of the distances between the nodes of the wave function. For the supersymmetric 
Hamiltonian it is not very convenient to consider the "Ricatti" variable z = — because it is not an 
additive process but a multiplicative process: dxZ = k + kz^ — 2(f){x)z. To help the discussion we 
introduce two intermediate variables, a phase variable 'd{x) and an envelope function exp^(x): 

^{x) = e«(^)sin??(x) (81) 
X(x) = -e«(^)cosi?(x) . (82) 

These two functions obey the set of coupled stochastic differential equations, written in the 
Stratonovich convention: 

— -d = k -\- (f){x) sva.2'd (Stratonovich) (83) 
dx 

= —(j){x) cos 2t? (Stratonovich) . (84) 

The fact that the periodicity in ■(? of these two equations is vr/2, and that V{x) is 5-correlated, 
means that the distance separating the two points where -d = nnr/2 and -d = (m + l)7r/2 depends 
only on V{x) between these two points (the periodicity in -d is 7r/2 and not vr because the sign 
can always be absorbed in the function (j){x) if it is a white noise of zero mean). Then the 
statistically independent random variables to be considered in a first step are not the distances i 
between successive nodes but rather the distances, denoted A, separating the points where '&{x) 
takes a value equal to a multiple integer of tt/2 (the point where ('!9(x)mod7r) = 7r/2 corresponds 
to a local extremum of the oscillatory part of the wave function). The length i between two 
consecutive nodes of the wave function is then the sum of two independent A's. 

It is more convenient to introduce a process which is additive in the noise (j){x). It is easy 
to see that this is achieved by the change of variable C{x) = ^ In | tan'i?(x)|. This variable obeys 
the stochastic differential equation: 

= ^cosh2C + (^(x) . (85) 

dx 

This is a Langevin equation for a "particle" of position C traveling from — oo to +oo in a potential 

i7(C) = -|sinh2C (86) 

(see figure |5|) and feeling the random "force" (/)(x). A is the "time" the "particle" needs to go 
through the interval. As for the model A, we are going to study the statistical properties of 
the distances i between the nodes, sum of two independent variables A. The knowledge of the 
distribution of i will allow us to find the distribution of the energies. 
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4.1 Distribution of distances between the nodes 



We introduce the moments of the "time" needed by the variable C{x) to reach +00 starting from 
the position 

KiC) = (A" I C(o) = C; C(A) = +00) . (87) 

They are given by (see appendix |A|): 



An(C) = -/ 

9 Jc 



(89) 



,w ^U{C') / ,^// --C/(C") A 
y j( J-00 

Ao(C) = 1. 

The moments of the random variable of interest are (A") = A„(— 00). 
4.1.1 Low energy limit: E <^ 

At low energy k <^ g, the potential is still a monotonic function and there is no process of 
trapping of the variable C by a well as for the low energy limit of model A, which means that 
we have to develop a specific approximation scheme. For this purpose we start by studying the 
average time Ai((^) to go from C to +00: 



Ai(C) 



2 

9JC 



+00 



T A/ - - sinh 2C' 

dC e 9 ^ 



dC 



(90) 



We introduce (^±, the two solutions of the equation exp ^|sinh2C^ = 0. In the low energy 

limit: C± — In f • The positions (± are the crossover points where the force deriving from the 
potential U{C) is of the same order as the random "force" (j){x). The study of (90) leads us to 
distinguish three regions to which the initial condition can belong to: 

(0 c+ < c 

(91) 



9 



(^^) C- < C < C+ 
[iii) C < C- 



Ai(C) 



(C+-C-)'-(C-C-)' 



+ 



Ai(C) ^ - ^ 
9 



1 



+ -(C+-C-)Vi. 

9 9 



(92) 
(93) 



3^ 



1 


v.. 








1 



Figure 5: Potential "seen" by the variable C, for fc = 10 
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What can we learn from these behaviours ? (i) if ^ starts from the neighbourhood of C+j it 
needs a "time" of order \/g to reach +00. (ii) When (" starts somewhere in the "time" 

needed to reach +00 is dominated by the first term, which behaves hke C^, i.e. hke the "time" 
for a free diffusive "particle", {iii) When C, starts from —00 it needs a "time" 1/g to reach 

2 



then ^ travels from C„ to in a "time" ^dTffuston 



~ C-)^ — I lii^ f and eventually ends 
at +00 after an additional "time" 1/g. 

The physical picture for the motion of the fictitious particle of position C, we get from these 
results is now very clear. Q travels from —00 to C- very quickly due to the potential only; for 
C £ [C- ) C+] the potential becomes negligible compared to the random force which is of order g 
and C, evolves due to the random force only. It increases like the position a free diffusive particle 
until it reaches C+ from where it goes very rapidly to +00. 

We are now going to use this picture to study the distribution of A. We are interested in the 
characteristic function for the traveling time 



Ma,C) = (e-"^IC(o) = C; C(A) 



-00) 



(94) 



From the Langevin equation (|8^ 
appendix |A|) 



we get the BFPE generator and we see that h obeys (see 
/i(a,C) = a/i(a,C) (95) 



with boundary conditions: 



5c/i(a,C)Loo 
/i(a, +00) 



0, 
1 . 



(96) 
(97) 



Since the traveling "time" is dominated by the "time" spent in the region , Q^] where the 
diffusion is free, it means that the characteristic function in the limit k <^ g is approximatively 
given by the solution of the equation for the free diffusion on this finite interval. Then equations 
(95, 96| , |97| ) may be replaced by the following equation 



with boundary conditions: 



dt;h{a,C-) 
/i(a,C+) 



a h{a, C) 



0, 
1 . 



(98) 

(99) 
(100) 



The first condition is a reflection condition since the variable C "sees" a steep wall in The 
second condition is an absorption condition since as C reaches C+, it eventually ends to +00 after 
a negligible time. Since we are dealing with free diffusion, the solution is then very easy to find: 



cosh 



h{a,C) 



KC-C-) 



cosh 



/2a /A 
-(C4 



The characteristic function for the length A is: 

(e""^) = h{a,C- 

where we have introduced B 
I between the nodes of the w 
give the characteristic function for t. 

-al\ 



cosh - 



(101) 



(102) 



In^^ 



r — TT- In^ . Since we are in fact interested in the distance 

g k 2g E 

i between the nodes of the wave function, which is given as the sum of two independent A's, we 



(e- 



cosh 



(103) 



16 



The inverse Laplace transform gives the distribution function: 



Pii) = I 



dq 



2i7r cosh yjq 



(104) 



where the integral is taken over a Bromwitch contour, cosh yfq is an analytic function in the 



variable g, whose zeros are: q„ 



-^{2m + 1)^ for m G N (note that there is no branch cut 



since cosh is an even function) . We can expand cosh in the neighbourhood of its zeros and 
we get: 



cosh ^ ~ 



(-ir 



g^qm 7r(2m + 1 



1 + 



1 



7r2(2m + 1)2 



(g - Qm) + 



(105) 



We use the residue's theorem to evaluate integral ( |104| ); we can indeed check that the contribu- 
tion of integrals over the large semi-circles needed to close the contour in the complex plane are 
vanishing in the limit of their radius going to infinity. Then we get: 



1 



Using the identities Q: Yl 



m=0 

1 



^_2L_(2m+l)2£/B 



m=0 (2m+l)2 



=0 (2m+l)4 ~ 96 

we can check the normalization and get the two first cumulants: 



III 1 
T — 96 '^^^^ Z^m=0 (2m+l)B 



(106) 



960' 



(107) 
(108) 



The average distance between consecutive nodes is the inverse of the integrated density of states. 
We can check from the exact result ( |7^ ) that the limit behaviour of N{E) is indeed in perfect 
agreement with the result we find here {£) = N{E)~^ ~ which gives us a certain 

confidence in the approximation we have made to calculate the characteristic function h{a,C)- 
Since B = {£) it can be replaced by N{E)^^ and we can rewrite the distribution as 



where 



P{£) = N{E)mo{N{E)e), 



wo{X) = e{X) [7r2(2m + ifX - 2] 

m=0 



'-{2m+lfX 



(109) 
(110) 



replaces the exponential function obtained for the low energy limit in model A 

We can extract the limiting behaviours of the distribution (|106D . For this purpose we intro- 
duce the ^-function (not to be confused with the Heaviside function 9{x)) 



e{x) = e- 



m=0 




1 + 2 ^(-l)'"e- 



4x 



an) 



m=l 



(this is the elliptic theta function i?i(7r/2 | Aix/tt) [g4|, 8.180]). The distribution ( pTOD is related 
to the 0-function by: 

2 12X^ + A9( —X^ . (112) 



ruoiX) 



dX 



Using (111, 112| ) we can find an expression adapted for the limit X — > 0: 



4 e{x) 



-i/x 



(113) 



17 



The tail of the distribution is exponential and dominated by the term m = in ( |110D : 

wo{X) ~ TT^Xe-^^ . (114) 




0.5 1 1.5 2 2.5 3 



X=l/B 



Figure 6: Function wo{X) giving the distribution of the distance between consecutive nodes of 
the wave function for the super symmetric Hamiltonian in the limit E ^ 0. Dashed: term m = 1 
of (|ri3|); Dotted: term m = of ([TlO|). 



Structure of the wave function at low energy 

It is well known that the supersymmetric Hamiltonian ( |76| ) exhibits a delocalization transition 
at zero energy [^2|. As energy goes to zero, the localization length behaves like — | ^T^i^d"^ /E). 
The distribution we have just found for the distances between the nodes shows that two con- 
secutive nodes are separated by a distance of order B = ^ \v?{g'^ /E) which is much larger than 
the localization length A. It is also worth mentioning that this distance is the correlation length 



appearing in the average Green's function |22]. By contrast to what happens for the model A 
at low energy, where the nodes of the wave function can be arbitrarily close, the nodes of the 
wave function for the supersymmetric model B with /i = are extremely unlikely to be closer 
than a distance of order | \r?{g'^ /E): the behaviour of the distribution P{i) (figure P) indicates 
a "repulsion" of the nodes of the wave function. 

4.1.2 High energy limit (£J ^ 9^)' small disorder expansion 

At high energy we perform the same perturbative analysis in the disordered strength g as for 
model A, hence we do not enter into the details. The generating function for the cumulants 
w{aX) = In /i(a,C) obeys: 

k cosh 2C dc_w + I [dlw + [dc^wf] = a . (115) 

We solve this equation perturbatively in g: w = w^^'^ +w^^'^ + • • where w*^""^ = 0{g"'). Since the 
reflection condition is at C = —oo and the absorbing one at (" = +oo, we have now: w{a, +oo) = 
to ensure h(a, +oo) = 1. To zeroth order we find: 

„<o.(„,c,=--/ dC'^ (116) 



and to order n: 



cosh2C' 



n — 1 



m=0 



(117) 
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The information of interest is contained in the first two orders: 

w{a,-oo) = -a^ + ^^ + 0{g^), (118) 

that give the cumulants of the variable A. Let us recall that the length £ is the sum of two 
independent A's, then: 

W = l + 0{g'], (119) 

Since the cumulants of higher orders are small, {{i^)) = 0{g^~^), the distribution of £ is 
Gaussian in the small disorder limit as for the high energy limit of model A, given by (^Op. 



4.2 Distribution of energy level 
4.2.1 Low energy: E <^ 

We first consider the distribution for the ground state energy and will give an explicit form 
for the excited state energy distribution without going further into the calculation to avoid 
technical increasing complexity. The probability for the ground state to be at a given energy is 
proportional to the probability that the distance between the two first nodes of the solution of 
the Schrodinger equation is equal to the length of the disordered region. According to (|6^ we 
have: 

Wo{E) = ^^P{i = L) . (121) 
Using (|109| ), we see that this distribution is a scaling function of the variable: 



X = L NiE) = — , (122) 

which is the averaged number of states below E for a system of length L. Then 

Wq{E) = Lp{E) vjo{L N{E)) . (123) 

From ( |ll3| , |Tl^) we see that the distribution Wq{E) presents a log- normal behaviour at low 
energy 

«'»<^'-^|-P-T^ E«.Vv^ (124) 

and the following behaviour at large E (however smaller than , not to be out of the range of 
validity of the approximation we have made): 

/e-^«^«.^ (125) 

(see figure |^ . 

We may now proceed to study several aspects of the distribution Wo{E). The typical value 
of the distribution is given by the limiting behaviour of the distribution for small E (124): 

- 5^e-3^ . (126) 
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Figure 7: Distribution Wo{E) for the supersymmetric potential with /x = 0. L = 10, g = I. 



It is also interesting to estimate the median value: Jq " dE Wo{E) = 1/2. Assuming that we 

'ln{g2ABg-d)^ 



can also consider the limit behaviour ( 124 ), we find that it is the solution of ^ 
7/8, where ^{z) is the error function. Then: 



(127) 



(c is a numerical factor) which is at the boundary of the domain where approximation ( [1241 ) 
holds. 

We may also compute the moments of the ground state energy: 

{E^)= dEE^Wo{E)= dXwoiX)e~^V'^3L/x _ ^^28) 
Jo Jo 

If /3 > the exponential select the tail ( |114| ) of the distribution: 

(E^) ~ TT^ r dX X Q-^^-f^^ ^ . (129) 
h 

This integral is easily worked out by the steepest descent method, and we eventually find: 

for ;3 > . (130) 



l^(5L)V2exp-^(vrW)^/^ 



~ e <^^^^^ which 



In particular, the mean value of the ground state energy behaves like {Eq 
is much larger than the median value. This result is in agreement with the work of Monthus 



et al. |3£] who found upper and lower bounds using a perturbative expression for the ground 
state energy as a functional of (t){x). The atypical behaviour of the positive moments with n like 
{Eq"") ~ i)"'^^^ may also be recovered with the method used in [^] as remarked by Oshanin p5[ . 

It is also possible to compute the negative moments. If /5 < 0, the exponential in (ll28| ) 
selects the origin of the distribution. The steepest descent method gives: 



L^oo 



exp -0^gL 



for /3 < 



(131) 



which presents a dependence in j3 characteristic of the log-normal behaviour (|124| ). 
We now make some remarks. 



(i) It is worth mentioning that the fluctuations are behaving at large L like 6Eq ~ exp — |(47r^5(L)^/'^ 
and are much larger than the mean value (Eq)'^ ~ exp — |(87r^(7L)^/^. Thus the ground state 
energy is not a self averaging quantity in the L — > oo limit as it was the case for the model A. 
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(a) A remark related to the previous one. Since the distribution Wo{E) is a scaUng function 
of L/ln^ (!/£;), it becomes more and more peaked near zero when L oo. Moreover it is 
rather obvious that it can not be written as a scahng function of (E — Ei{L))/E2{L) where 
Ei^2{L) would take into account the L-dependence as it was possible for model A. In other 
words, none of the moments ( |130| ) determines the nature of the fluctuations at L ^ cxd since 
{Eq)^/^ ~ L^/^'^ exp — |(7r^L//3)^/^ has a non trivial dependence in j3. 



[in) It is interesting to note that the average Green's function, derived in [22|, presents the same 
kind of behaviour as the distribution Wq{E): 

{x\^^jYp^\x')j ~ EmCmexp (-^^(2m + l)2|a; -x'M (in the limit ^ ^ with g = 1). 



{iv) The distribution ( |11[1| ) was obtained in |^3[ in the context of the classical diffusion by using 
a real space renormalization group method. In this case the distribution is interpreted as the 
distribution of the smallest relaxation time. 



[v) The distribution ( |110| ) has still another interpretation: it is related to the distribution of the 
span of a Brownian motion (see ||3^ and references therein). 



{vi) The n-th excited state energy distribution is related to the distribution 7^„(A) of the sum of 
the n + 1 lengths I by ( |68| ) and may be studied by the same kind of calculations, starting from: 

/■+ioo J ^qLN{E) 



where we have used the fact in (103) that B = {£) = N{E) ^ = In^ (g"^ / E) . For example we 



find that the first excited state energy distribution, Wi{E) = Lp{E) wi{L N{E)), involves the 
scaling function: 



m=0 



— {2m + 1)^X3 - 2^2(2m + 1)^x2 
6 



(133) 



eix) 



oo 



where X is given by ( |122| ). Its limiting behaviours are: 

32e{X)_^-4/x 



J2 {-irm\m^ - l)e-"^ , (134) 

m=l 



(135) 



vr 



4 



,2 



_ —X-^e--^ . (136) 

X^oo 6 

Note however that for large n, since the moments of £ are finite, we expect a Gaussian 
distribution for Vni-^) due to the central limit theorem. 

(mi) The distribution Wq^E) is not the distribution of the extreme value of a set of independent 
random variables (see appendix]^) as it is the case for model A. We will come back to this point 
at the end. 

4.2.2 High energy: E ^ 

For high energy the analysis is very similar to the one that was done for the model A since the 
distribution of lengths ts is the same. Only the behaviours of several quantities change. The 
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distribution Wn{E) is given by (|70| ) with s{k) = (n+ The condition to be at high energy 

leads to the same expression (^Tj) with the only difference that the localization length for the 
supersymmetric model at high energy reaches a constant value: A ~ |. Then ( [7^ ) still holds 
with the mean value of the energy still being the free result (|73|) and the fluctuation being 



5En = {n + l)^J^. (137) 

The fluctuations now depend on n, compared to model A, and its L-dependence also changes. 
However they remain small compared to the mean value (En) due to condition (|7l|). It is also 
interesting to note that the ratio of the fluctuation and the mean level spacing has the same 
form: 

6El 2 L 



which means that the absence of level repulsion may only be expected for the localized regime. 



4.3 The case /x 7^ at low energy (E g^) 

We consider in this last section the case where the function (l){x) is distributed according to ([77|) 
with a non zero mean: {4>{x)) = fig (note that /i is a dimensionless parameter). We will only 
discuss the low energy case since we are not expecting any modification at high energy. We have 
to come back for a while to equation (^) for the phase from which we want to get an equation 
for an additive process. Performing the change of variable (" = ^ln|tani?|, we have in fact to 
distinguish two cases: if (?9 mod vr) G [0,7r/2] then l/sin2'!9 = +cosh2(^ and we get (^), but 
if ("i? mod vr) £ [7r/2,7r] then 1/ sin 2?? = — cosh2C and we arrive at d^C = — A;cosh2(^ + (j){x). 
The change of the sign in the potential is related to the fact that C, decreases if ^ increases if 
(-i? mod vr) E [7r/2, vr]. So a more convenient change of variable is in fact: 

C = ^ln|tani?| if (1? mod vr) G [0, 7r/2] (139) 

C = -^ln|tanT9| if (i9 mod vr) G [7r/2, vr] . (140) 
With this convention, is always traveling from —00 to +00. Then obeys 

= fecosh2C + (/'(x) if (t? mod vr) G [0,7r/2] (141) 

ax 

and 

= k cosh 2C - 0(x) if {'d mod vr) G [tt/2, vr] . (142) 

ax 

It follows that we have to consider alternatively two different equations. We did not mention 
this point before because if (p^x) is a white noise of zero mean, since ( decorrelates when it 
reaches +00, the difference in the sign can always be absorbed in the white noise. We rewrite 



these two equations in terms of a normalized white noise rj{x) of zero mean: (??(x)) = and 
{r]{x) r]{x')) = S{x — x'): 

= A;cosh2C + /UC/ + ^r?(a;) if (i9 mod vr) G [0, 7r/2] (143) 

= k cosh2C - fig + ^/gr]{x) if (i9 mod vr) G [7r/2, vr] . (144) 
We have to introduce two potentials ?7i(C) and U2{C)- 

UiiC) = -/i<7C-|sinh2C (145) 

^^2(0 = +wC-|sinh2C, (146) 
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plotted in figure ra. 
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Figure 8: The two potentials in which the variable C evolves alternatively. 

We must also introduce two variables Ai and A2 that give the "times" the variable C, needs 
to go from —00 to +00 in potential Ui{C) and f/2(C)) respectively. For a drift fj, = those two 
variables have the same statistical properties as it was implicit in the previous section but if 
// 7^ 0, their distributions are different. We are now going to study these two distributions, 
denoted pi(Ai) and P2{^2)- We introduce the two characteristic functions for the "times" Ai^2 
needed by the random process obeying ( |143|J144| ) to reach +cxd starting from C: 



/ii,2(a,C) = (e" 



-aAi 



C(o) = C; C(Ai,2) = +00) 



(147) 



We recall that /ii,2(a, —00) = dApi^2(A) e These two functions obey (see appendix 



{k cosh 2C ± fj,g) 5^ + ] /ii,2(a, C) = " ^i,2(a, C) 



(148) 



where the upper sign (here +) corresponds to 1 and the lower sign (here — ) corresponds to 2. 
The boundary conditions are as usual d(hi^2io(, C)\-oo — ^1,2(0, +00) = 1. Since we are 
dealing with the limit — > 0, we can make the same approximation as for the /j, = case, i.e. 
replace the previous equations by an equation for the free diffusion on the interval since 
the time spent out of this interval is negligible. Then we have to solve: 



^±^ic/5^ + l^l] /ii,2(a,C) = a ^1,2(0,0 
with the boundary conditions: 



^1,2(0, C+) 



0, 
1 . 



The solution is easily found: 
^1,2(0,0 = 

where 



.TM(C~a) - C-) ± ^ smh7(C - C-) 

cosh 7(C+ - C- ) ± ^ sinh 7(C+ - C- 



7 



2a 



The characteristic functions for Ai^2 are: 
-aAi 2 



(e- 



^1,2(0, C- 



g±M{C+-C-) 



cosh 7 ((^4 



C-)±^sinh7(C4 



(149) 



(150) 
(151) 



(152) 



(153) 



(154) 
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We can more conveniently consider the generating function tt;i^2(cK) C-) = li^ ^1,2(0, C~) for 
the cumulants of Ai^2- In the hmit of low energy, when the condition 



^(C+-C-)=/"ln5/fc»l (155) 
is fulfilled, a careful analysis shows that: 

w\\a, C-) = -a h — ^-5 VO \a (156) 



2! p?g^ \ ^^g^ 

and 



-2(a,C-) = - In (1 + ^e^MC^-^-)) + ^ k-"'" ^^^L^"^ • (1^7) 



The characteristic function ( |156D corresponds to a sharp Gaussian distribution p\ ( Ai ) for 

(Ai) = (158) 

((a;» = (159) 

This result is not surprising since the potential \J\ is monotonic and linear on the interesting 
interval (see figure |^, we expect the "time" to go through the interval to be '^gp^^'^ = ^^~^g~ ^nd 
to fluctuate weakly for large enough drift. The precise criterium for the validity of this result is 
that the time characterising the motion due to the drift, Tdrift = is much shorter than the 

time characterising the motion due to the diffusion, Tdig = ^d'jffuskm ' ^'^^^ ^'^^ i^a^ , 
which leads to (11551). 



In (157), the first term gives contributions to the cumulants exponentially large in ((^+ — QJ) 
whereas the sum gives contributions powers of (i^+ — C_). Then we can forget the second 
term and the characteristic function is the logarithm of the Laplace transform of a Poisson law 
P2(A2) = ■^exp-A2/(A2) with 

The fact that the distribution P2(A2) is a Poisson law was expected from the shape of potential 
[/2(C) (see figure 1^: the potential possesses a well that is able to trap the variable C for a long 
time. We can now compare the result of the approximation we have made for (£) = (Ai + A2) 

with the exact result (|7|). We have found [1) ~ ^^^^^ ^ ^g^k)^^ - ^g{^^ whereas (||) 



gives N{E) ^ ~ ^ sin^ n^i ' I-'^^Pi^^ ^he pre-factors being different, our approximation is 

indeed able to give the well-known power law behaviour of the integrated density of states |25, 22 1: 
N{E) ^ E^^. We may also have used the Arrhenius formula (|193|) to find (A2) ~ — ' ^1 



fj-g \ k ) ' 

which presents still a different pre-factor (equation (|193| ) does not give the correct pre-factor 
maybe because the potential ^^2(0 is not smooth enough in the neighbourhood of its local 
minimum in the limit A; ^ as it is assumed to derive the Arrhenius law ( |193| )). However we 
may distinguish two levels of approximation in what we have done, (a) We have shown that the 
distribution of A2 is Poisson, that could be demonstrated in a more general way following the 
proof presented in appendix ^ for the case of a potential possessing a local minimum, {h) We 
have given an approximative expression of the average time (A2) that only gives the correct 
behaviour with E but not the correct pre-factor. 
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However we can avoid the not absolutely satisfactory approximation (b) because we know 
that (Ai + A2) = N{E)~^ is exact. Then we may give the distribution of the length £ between 
nodes of the wave function: 

Pi£) = [ dA2P2(A2)pi(^-A2); (161) 



we have shown that P2{^2) varies on a characteristic scale (A2) ^ (Ai) ^ ii^i)) compared 
to which pi(Ai) is very narrow. It follows that P{£) — P2i£) apart from the behaviour at the 
origin we are not interested in, since it is associated with small samples. Then we conclude that 

P(^) = iV(^)e-^^(^) . (162) 

We can now give the distribution of the n-th excited state which has the same form than 
what was found for the model A (H9). Then we have: 



Wn{E) = Lf,a,E'^-' ^^''^p\ -^'^^^' , (163) 

n! 



where the coefficient is defined by: N{E) ~ a^E^. We recall that this result is valid if ( |155| ) 
is fulfilled, that is filn(g/k) ^ 1, in addition to the fact that k <^ g. 
It follows that the distribution involves the scaling function: 

WniE) = 1 u;„ (^) (164) 

where 



Se \ 5e 



vn/i 

a;„(X)=^X'^-i^e-^', (165) 



n 



the energy scale being 



5e = ^-^ . (166) 

We may easily compute the moments: 

T[n + I) 

It is also interesting to compare the mean level spacing (A„) = {En+i — En) and the fluctu- 



ations 6 En = yJ{{El)). We get: 



m/^x + n + l) 

A„) = 6e — — — . 168) 

^l{n + 1)! 



Then 



I , r(2/^ + n + l) 



(A„) r(l/^ + n + l) 

which becomes large at small ^ (^ l/v?): 



(A„> " 

At large /i the ratio reaches a constant value: 

SEn 



n/2+5/42l/At _ {170) 



hm — ^ = (n + l)VV'(n + l) , (171) 

^l^oo {An) 
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where ip{z) is the digamma function (in particular = 7r^/6). For small n the fluctuations 

and the mean level spacing are of same order. 



The distribution (163) has the form of the distribution of the n + 1-th lowest variable among 
— > oo statistically independent variables distributed by a law behaving like p{E) ~ E^^^ (see 
appendix]^), the same behaviour as the density of states. This remark shows that the energies 
behave like statistically independent ordered variables, in agreement with the expected absence 
of level repulsion in the localized regime. 

5 Conclusion 

We have considered two one-dimensional disordered models: one with diagonal disorder (A) and 
another with off-diagonal disorder (B). 

For model A we have derived the distribution of the distance between consecutive nodes of 
the wave function in the limit \E\ ^ o"^/^ both in the negative part of the spectrum and in the 
positive part. Using these results we were able to find the distribution of the n-th excited state 
(|4^|7^ ) . For < we have shown that ( ^9|) is a scaling law (|6ll ) , similar to the extreme statistics 
of independent variables. If E < the typical value of the n-th energy behaves with the size of 
the system like: E^"^ ~ —In^^^L. The width of its distribution behaves like 6En ~ In^^^'^L. If 
-E > we have found E!^'^ ~ 1 /L^ and 6En ~ 1/ Note however that the relative fluctuations 
of En in this latter case are small despite the behaviours with L suggest the opposite. 

For model B we have first considered the case /U = (mean value of the function entering 
the supersymmetric potential). The high energy limit (universal regime) gives results similar 
to those for model A. In the low energy limit we have found the distribution for the ground 
state (123, llCl| ). We have shown that this distribution is broad, its positive moments being all 



dominated by the tail of the distribution: E'J" ~ e"^ <C E;^"'^ ~ e"^ < (£^o> ~ e-^^^\ 
The moments have an atypical n-dependence: (Eq) ~ 6"^''^. We have also given explicitely the 
distribution for the second energy level ( 133|jl34| ) and an integral representation for the other 



energy levels (132). For ^ = these distributions do not have the form of the distribution of 
extremes of independent variables. For ^ 7^ we were able to derive the distributions Wn{E) in 
the low energy limit (in the high energy limit, universal regime, we do not expect any difference 
with the picture obtained for fi = 0). We have shown that Wn{E) is a scaling function oJn{E/5e) 
where the energy scale behaves like: 6e oc L~^/^. For /i 7^ the distribution exhibits extreme 
value statistics as for model A. 

We now discuss the relation between the distributions Wn{E) we have found and the extreme 
value statistics (see appendix |^. We have seen that for model A and for model B with ji 7^ 
0, Wn{E) has the form of a distribution of extreme values. If we suppose that energies are 
behaving like statistically independent random variables, the distribution of one of them should 
be proportional to the density of states; the number N of these variables is proportional to 
the size of the system. If we replace by L x and p{x) by p{E)/^} in ( |196D , where $7 is the 
total number of states per unit length (in principle infinite since the spectrum is unbounded from 
above), it is easy to see that we get the equation ( ^9|) obtained both for model A and for model B 
with /U 7^ 0. We now consider more specifically the case of the ground state (n = 0), since Wq(E) 
is directly related to the distribution P{i), and proceed in the opposite way we have followed until 
now. If we admit that the distribution of the lowest energy is Wq{E) = Lp[E) exp —L N{E) (see 
equation (|49D), as a consequence of the statistical independence of the energies, then we conclude 
that P{t} is Poissonian. Conversely if P{1) is not Poisson in the range of energy where Eq is 
expected to be found, there might exist some correlations between the energies. This is indeed 
the case for the high energy limit: the distribution ( [72|) for n = 0, a consequence of the narrow 



Gaussian distribution P {(.)■, is valid only in the delocalized regime due to condition (71), and in 
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this regime level repulsion occm's, as explained, due to ([751) . This argument can only be used 
for the ground state energy distribution since Wn{E) is not directly related to P{i) in this case. 
So the fact that ( [72[ ) for n ^ 1 comes from a Gaussian narrow distribution P{i) in this range 
of energy does not necessarily mean level repulsion. Since the distribution of the ground state 
for the super symmetric model with /i = is Wo{E) = Lp{E)wo{L N{E)) with wo{X) ^ e""'^, 
this suggests that level repulsion might occur at the bottom of the spectrum. This seems to 
contradict the expected absence of level repulsion in localized regime since Wq{Ei) shows that 
the most probable energy E ^ ~ exp —gL, i.e. A ~ ^ln{g'^/E) <C L, are associated 
with localized states. However A ~ ^ln{g'^/E) is the localization length for the infinite size 
system, given by the inverse of the mean Lyapunov exponent (7) = 1/A which is going to zero 
if E ^ 0. For finite size system the Lyapunov exponent has some Gaussian fluctuations which 
cause fluctuations of the localization length: ((7^)) oc 1/L. For a given L the fluctuations of 
the Lyapunov exponent become of the order of its mean value when the energy becomes smaller 
than an energy g^ exp(— c-y/ gL) of the order of the median value of the distribution Wo{E); in 
this case, the large fluctuations of the Lyapunov exponent may cause delocalization and it might 
not be surprising that level correlation appears in that regime being the reason of the fact that 
Wo{E) does not behave like the extreme value distribution of independent variables as it is the 
case for /i 7^ or for model A. Nevertheless a deeper understanding of level correlations seems 
to be needed for the supersymmetric case with fi = 0. 

Since model B is related to the problem of classical diffusion in a random medium, it would 
be interesting to know if these results have an application in this case. The energies of the 
Hamiltonian should be related to relaxation times [^]. The spectral properties near E = are 
important for the disordered spin chain models |21], the distribution Wo{E) of the lowest mode 



of the Hamiltonian (|7q ) might also have some interest in that context. 
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A Distribution of the escape time for a diffusive particle trapped 
in a well 

In this appendix we recall known results about the trapping of a Brownian particle by a well 
and the distribution of the escape time. In the regime of interest here the average escape time 
is given by the Arrhenius law. Most of what is in this appendix may be found in standard 
textbooks like l47| ] and we summarize here some ideas needed throughout this article. 

We consider a particle whose position x{t) obeys the following stochastic differential equation: 

dx{t) = A{x)dt + y^B{^dW{t) (Ito) (172) 

where W{t) is a normalized Wiener process: {W{t)) = and {W{t)W{t')) = min(t,t'). This 
equation is understood in the Ito sense. The propagator p{x', t\x, 0) for the diffusion (conditional 
probability for the particle to be at x' at time t, starting from x at initial time 0) obeys the 
backward Fokker-Planck equation (BFPE): 

dtp{x', t\x, 0) = p{x\ t\x, 0) , (173) 
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where 

Gx = A{x)dx + \B{x)dl (174) 

is the BFPE generator. 

Let us consider a process starting from a point x(0) = x belonging to the interval [a, h]. Let 
us call T the time needed by the particle to leave this interval. The probability that the particle 
is still in the interval after time T is 

/•b j-oo 

/ dx'p(x',r|x,o) = / dr'p^(r'), (i75) 

J a Jt 

where Px{T) is the probability density for the escape time. Let us define Tn{x) = (T" | x(0) = 
x; x{T) = a or b), the n-th moment of the escape time from the interval [a,b]. It follows from 
(|T75|) that 

GxTn{x) = -nTn-i{x) . (176) 

Since the particle may escape from both sides of the interval, we have to impose the following 
boundary conditions: Tn{a) = Tn{b) = 0. This absorbing boundary conditions mean that the 
particle leaves the interval and never comes back as soon as it reaches one of the two edges. In 
this sense we are interested in a time of first exit, relevant for the question considered in this 
article. 

In the following we will be interested in the more simple case of a reflection condition at one 
side of the interval, x = a, so that the particle may only escape from the side x = b: 

Tn{x) = (T" I x(0) = X- x{T) = b) . (177) 



The reflecting boundary condition for the BFPE is dxp{x' ,t\x,Q)\x=a = ||4^ which implies the 
following boundary conditions for the moments: 

dxTnia) = (178) 
Tn{b) = 0. (179) 

In this case it is easy to construct from ( 1761) the moments: 
where 

r ,2A{x') , , 

V'o(a;) = exp / dx „, . (181) 



B{x') 

Keeping in mind that Tq{x) = 1, it follows that the moments may be computed recursively. 
We also introduce the generating function for the distribution Px{T\. 



h{a, x) = (e""^ I x(0) = x, x{T) = b) , (182) 



which may be used to analyze the distribution. It is clear from ( |176| ) that it obeys the following 
differential equation: 

Gxh{a^x) = ah[a,x) . (183) 
The boundary conditions for reflection in a and escape (absorption) at b are: 

dxh{a,a) = (184) 
h{a,b) = 1. (185) 
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Instead of considering (|183| , |184|Jl8^ ) it can be more convenient to write an integral equation for 
h: 

1 r' , „ Mx") 



h{a,x) = 1 — 2a / dx' 



dx" 



h{a, x") . 



(186) 



V^o(xOya ""^ B{X") 

We now derive the distribution Px{T) when the particle is trapped by the well of a potential 
U (x), in the small noise limit. We consider the more simple situation of a stochastic differential 
equation of the form: 

dx{t) = -U'{x)dt + Vd dW{t) (187) 

where D is the diffusion constant. The shape of the potential of interest is depicted in figure 
|9[ The potential U (x) has a local minimum xq from the neighbourhood of which the particle 
starts. Let us call Xi the initial position of the particle. The potential has at xi > xo a local 
maximum. We are interested in the time that the particle need to escape the well, jumping over 
the potential barrier due to a fluctuation. The interval [a, b] to be considered has to include the 
well and the barrier. As it will be clear in the following we have to consider a limit b such that 
the distance b — xi is sufficiently large (the relevant length scale is given by the curvature of 
U{x) at xi). For similar reasons the distance xq — a has to be sufficiently large as well. Apart 
from this two restrictions the precise positions of a and b are of little importance. 




Figure 9: The particle starts from Xi and escape the well after a time T. 

Let us analyze the behaviour of the moments when the diffusion D is small compared to the 
height of the barrier. Equation ( |180| ) now takes the form: 

(188) 



T„(x,) = ^ Tdxe^^^^)/^ rdx'e-2^(-')/^r„_i(x') 



From the shape of the potential of figure ^, it is clear that the integral over x is dominated by 
the neighbourhood of xi and the integral over x' by the neighbourhood of xq. Accordingly, we 
have: 



r„(x,) ^r„„i(xo) / dxe2^(-)/^ / ' dx'e-2^(-')/^ . (189) 

Since we are dealing with the limit D — > 0, the two integrals may be estimated by the steepest 
descent method. Expanding the potential in the neighbourhood of its two extrema: 

U{x) C/(xo) + ^^^|^ (190) 

x^xo /Oq 

U{x) t/(xi)-fc^, (191) 

x^xi ZO^ 



b rxi 



we eventually find: 



r„(xi) ~nr„_i(xo)27r 50(^1 ^ . (192) 
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For n = 1, we recover the well-known Arrhenius law |47]: 



„ U(x-, )-U(xn) 

ri(xo) ~ 2^(5o5ie2 5 . (193) 

Let us remark that for 6 = xi we find half of this result. It is now straightforward to see that 
the moments are: 

Tjxo)~n! [ri(xo)r , (194) 

i.e. the moments of a Poisson law. It follows that the escape time from a deep well is distributed 
by: 

P,.,(r) = -f-e-^/^^(-). (195) 

Let us stress that the initial condition Xi played no role in the previous discussion: Tn{xi) ~ 
Tn{xo)- 

The same analysis may also be performed using (|186| ). 

B Extreme value statistics 

We give here a brief discussion on extreme order statistics, details of which can be found in 



|41|, g7 



We consider a set of M statistically independent variables Xn, distributed according to the 
law p{x). Then the M variables are ordered and we call Wn{x) the probability density of the 
n-th variable among these variables. We have: 

(px \ n— 1 / f+oo \Af—n 

/ dx' p{x') I p{x)dx I / dx' p{x') I . (196) 

We are now interested in the behaviour of this expression in the limit M ^ oo. Three kinds 
of distribution are usually distinguished. 

Type I. p{x) is unbounded from below and decays exponentially. In the limit N ^ oo the 
distribution Wn{x) is expected to be peaked around a value — oo so that we may use 

the fact that p{x) may be locally approximated by an exponential in the neighbourhood of x^^. 
Then we get, up to a rescaling Wn{x)dx = Vn{y)dy: 

Vn{y) = -i -TT exp (ny - ne^) , (197) 

where y is the rescaled variable. The relation between x and the scaling variable y depends on 
the distribution p{x) but not the scaling function Vn{y)- 

Type II. p{x) has a power law tail. We do not discuss this case which is not relevant for what is 
done in this article. 

Type III. p{x) is bounded from below. We assume its support starts at x = and it behaves 
like p{x) ~ x^~^ for small x. Then, up to a rescaling of the variable, we end with: 

vn{y) = -J^y^^-^e-^y' . (198) 
n — 1 ! 
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